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Abstract 

A dynamical theory which incorporates the electron-electron correlations and the 
effects of external magnetic fields for an electron escaping from a helium surface is 
presented. The degrees of freedom in the calculation of the escape rate is reduced from 
3A^ to 3 as compared with other approach. Explicit expressions for the escape rate 
in various situations are obtained. In particular, in the weak parallel magnetic field 
limit the tunneling rate has an exponential dependence quadratic with magnetic field 
strength and an unusual exponential increase linear with temperature. 

PACS#s: 73.40.Gk; 71.45.Gm; 73.20.Dx 
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When many electrons sit in a metastable well near a helium surface, the escape of an 
electron from the well is no longer a single particle problem due to the Coulomb interaction. 
Electron-electron (e-e) correlations play an essential role in the escape process and the physics 
becomes very rich. This system is ideal to test our understandings of the escape from a 
metastable well when many body effects are important. A number of experiments have 
been performed on the escaping of electrons from a helium surface.?? Recent theoretical 
studies have mainly concentrated on how the static e-e correlations affect the escape rate. 
A comprehensive understanding has been obtained in treating the e-e correlations as either 
instantaneous following the motion of an escaping electron or not following it at all.?? The 
real dynamical nature of the e-e correlations, however, has not been explored adequately. 
As one can see from a rather similar problem of the escaping of a particle from a metastable 
well in the presence of an environment, the dynamical response can be dominant.?? In the 
present paper we develop a dynamical theory for the escaping of an electron from a helium 
surface to account for the effects of both static and dynamical e-e correlations and magnetic 
fields. The original 3A^ degree freedom problem?? is simplified to a 3 degree freedom one, 
with N the total number of electrons in the problem. Expressions in various limits for the 
escape rate in the terms of temperature, 2-d electron density, and external magnetic fields 
are obtained. Particularly, an unusual exponential increase linear with temperature of the 
tunneling rate in the presence of a weak parallel magnetic field is found. The dependence 
on the magnetic field is quadratic in exponent in this limit. 

We consider the experimental relevant situation in which the lifetime of the metastable 
state of an electron is much longer than the relaxation time of the 2-d electrons and the 
density of 2-d electrons is low such that the Fermi temperature is the smallest energy in 
the problem. The escaping events are then statistically independent of each other and the 
exchange effect of an escaping electron with 2-d electrons can be ignored. A separation 
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between the escaping electron and the remaining 2-d electrons for each escape event can 
be made. For simplicity, we shall ignore the weaker interactions of the escaping electron to 
the surface waves of liquid helium and the helium vapor atoms, which has been discussed 
elsewhere??. 

We start with the derivation of the effective Hamiltonian for the escaping electron classi- 
cally. This Hamiltonian will be the base to study the escape process. For quantum tunneling, 
an imaginary time path integral method following Ref.[13] will be used. The classical equa- 
tion of motion for the escaping electron is: 

m^Kt = -VVoizt) + X Be. + eE,„(Ri) , (1) 

aV c at 

with Vo{z) = Vw{z) + Vi{z) + Vn{z). Here m is the mass of an electron, e is the negative 

electron charge, c is speed of light, and is the hardwall potential, Vw — oo ior z < 0, 

= for 2; > 0. This potential mimics the fact that it costs energy ~ leV for an electron 

to go into the liquid helium, which is very large in the present problem. The image potential 

Vi due to the polarization of helium liquid is Vi{z) = —e^A/z, with A = (e — l)/4(e + 1). The 

dielectric constant of liquid helium e — 1.057. The potential is the total electric potential 

produced by the external applied electric field (perpendicular only) and the electric field 

produced by the mean density uq of the 2-d electrons, 

Vn{z) = -e[E,^ + 27re(l - AK)no]z . (2) 

The condition for 2-d electron escaping to z — 00 from the surface is £^ea;+27re(l — 4A)no < 0. 
The external applied magnetic field is Bex, and Ei„ is the electric field produced by the 2-d 
electron density deviated from uq. 

The induced electric field Ej„ is generated by the 2-d density deviation from the mean 
value no, which in turn is induced by the escaping electron. Consequently, Ej„ can be 
expressed in terms of the motion of the escaping electron. The procedure is as follows. 
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Let the 2-d electron fluid sit on the surface of the liquid hehum, the x — y plane. The 
fluid is described by a set of hydrodynamical equations: the continuity equation and the 
Euler's equation. In the small density deviation and the nonrelativistic limit, we linearize 
the hydrodynamical equations. We solve for the density deviation, which is presented by 
the motion of the plasma modes. Then using the Poisson equation we obtain the induced 
electric field as?? 

/f tlqG'^ 1 f dt' 

J ^ u^^-u^Uk) J ^ ^-P{-^^*'} ^-P{-^(k . r, - u;t')} 

(1 - 4A)2 ex.p{-kzt} exp{i(k • - cot)}. (3) 

Here the plasma dispersion relation ujp{k) in eq.(3) is 

ujl{k) = a;| + ^27rA;(l - 4A) + (4) 
m m 

and the cyclotron frequency ujb = ei?_Lea;/mc with B±ex the component of the external 
magnetic field perpendicular to the helium surface. In the calculation the pressure p — nksT 
for the 2-d classical electron fiuid phase has been used. 

The induced electric field is contributed by the response of the environment, the plasma 
modes of 2-d electrons, to the motion of the escaping electron. Hence we are dealing a 
problem similar to the one in the discussion of the macroscopic quantum effect??, where 
the total Hamiltonian has three parts, a dissipative environment consisting of harmonic 
oscillators, a system of interesting, and the coupling between the system and the environment. 
Using this analogy, we find that the following effective Hamiltonian is equivalent to eqs.(l,3) 
to describe the motion of the escaping electron: 

1 o., . 1 



2m 



ex 

c 



+ Va(z)+ f dkY^ 

Jk<g^/no j=l,2 



^p^.{k) + -mu;l{k)x 



2 

^ /, N e^/n^(l — 4:A) ex.-p{—kz}cj\ 

0^) 27r\ 

^ mwp{k) J 



(5) 



with 

Va{z) = Vo{z) - ]- r^dk 27rfc^^5^%^i^ exp{-2kz} , (6) 
2 Jo mwp{k) 

and Ci — cos(k-r), C2 = sin(k-r), Bg^ = V x Aea;. Here gisa. numerical factor of order unit to 
be discussed below. The choice of the effective Hamiltonian is chosen in a form as in Ref . [14] 
in the discussion of the dissipative bath in quantum tunneling. When 2-d electrons follow 
the motion of the escaping electron completely, the so called the adiabatic limit, the response 
of 2-d electrons is described by the adiabatic potential given by eq.(6). The deviation from 
the adiabatic response is described by the plasma dynamics, the last term in eq.(5). The 
second term of eq.(6) corresponding to the correlation potential discussed in Ref.[l]. Hence 
an alternative justification of its usage in Refs.[l-5] is obtained here. Several features of the 
effective Hamiltonian should be pointed out. The coupling between the escaping electron 
and the plasma is highly nonlinear in the coordinate of the escaping electron. The damping 
of the escaping electron due to plasma is clearly superohmic case, when B±ex 7^ 0.?? It is 
also superohmic when B±ex = 0, a situation similar to the polaron problem. The effective 
Hamiltonian contains a weak temperature dependent through the plasma frequency, which 
comes from the equilibrium state of the electron fluid. The potential Va{z) is influenced by 
B±ex through the plasma frequency cub dependence. 

Because the average distance between electrons is l/^/no, there is no plasma mode of 
large wavenumbers k » ^/uq. Furthermore, the 2-d electron density deviation is large in 
this regime and the linearization approximation leading to eqs.(3,4) is not accurate. Then we 
need to introduce a cutoff plasma frequency ijJc = sj^ixgn^'^ jm corresponding to /c ~ g^fn^ 
which determines the fastest response of the plasma, as shown in Ref. [16] . This would suggest 
that we could only determine the numerical factor g in eq.(6) to be an order of unit, and 
the radius of the hole created by the escaping electron pushing 2-d electrons sideward?? 
to be the order of I/^/tiq. However, g can be theoretically determined accurately in the 
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following way: Let the escaping electron sit in the center of the hole created by itself in 2-d 
electrons. Because of the rotational symmetry there is no net force from 2-d electrons acting 
on the escaping electron. The escaping electron only feels the externally apphed force, i.e., 
—'^VAlhoie center = ^E^xZ. The choicc of g must satisfy this condition. Therefore we have 
g = 2y^. This completes the derivation of the effective Hamiltonian eq.(5). 

The effective Hamiltonian may have a wider application region than the region of validity 
of the hydrodynamical approach, so long as elementary excitations, such as the plasma modes 
calculated above, dominate the response of the 2-d electrons to the escaping electron. For 
example, this may include the case of the Wigner lattice phase. In this case, the plasma 
modes will be replaced by the phonon modes. One can expect that the form of the effective 
hamiltonian is the same as eq.(5) in the long wave length limit when B±gx — 0, because 
the density deviation is the same longitudinal one. As the local structure of 2-d electrons is 
not changed in the Kosterlitz-Thouless transition, we then expect the same short distance 
behavior. Then there should be no change of the adiabatic potential therefore no change of 
the escape rate cross the melting temperature. This is consistent with a recent experiment.?? 

We now calculate the escape rate starting from eq.(5) for various situations. In the high 
temperature regime, the escape is dominated by the thermal activation. The escape rate 
is = A exp{—£i,/kBT}. Here the prefactor A is weakly temperature dependent and is 
proportional to the density of 2-d electrons. The barrier height Sb is determined by the 
equation Sh — VA{zmax) ~ ^o, where z^ax is the position of the barrier top. The metastable 
(ground) state energy £o is directly calculated from the effective Hamiltonian of eq.(5) as 



electrons have been ignored. If we ignore its small effect on the ground state energy, there 





with q'q — Ti? /me 



:^A. Much weaker influences from the Stark shift and vertical spread of 2-d 
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is no influence of the parallel magnetic field S||ea; on the thermal activation rate. On the 
contrast, the activation rate will be increased in the presence of a perpendicular magnetic 
field B±ex, because B±ex affects the potential Va{z) in a special way. Specifically, we find 
that the barrier height Si, in the weak field limit as: 

^bW..=o - (^) In (^) (1 - 4A)2e2 , w » 1 • 

This result may be interpreted as that B^^^x effectively reduces the interaction between the 
2-d electrons and the escaping electron. For a strong magnetic field we should pointed out 
that the short magnetic length scale introduced by the strong magnetic field may indicate 
the invalid of the present hydrodynamics approach. 

In the low temperature region quantum tunneling dominates the escape. The tunneling 
rate is Tq — Ua exp{— 5'c/^}. The prefactor Ua is of an order of Eq/U. The classical action 
Sc = 5'e//[Rc(T)] is evaluated at the classical trajectory Rc('r) determined by the equation 
(55'e//[R(T)] = which is a 3 instead of 2A'" + 1 dimensional partial differential equation, with 
the effective action Seff as?? 

5e//[R(r)] = 5o[R(r)] + ^/dr / dr' J dk ^°^'j^^^~ exp{-u;p(^)|r - r'|} 

{[e^p{-kz{r)} - e^p{-kz{r')}f + 2 exp{-A;(^(r) + z{r'))}[l - cos[k • {r{r) - r(r'))]]} , 

(9) 

with 



So[R{r)] = Jdr 



^mR2(r) + i-Ae, • R + Va{z{t)) 
2 c 



(10) 



The nonlocal term in time in eq.(9) is a result of the reduction of the degrees of freedom 
from 3N to 3. Numerical calculation can be used in order to make a detailed comparison 
to experiments. In the following we discuss some prominent features of the tunneling rate 
analytically. 



First, we consider the case of the zero parallel magnetic field B\\ex — 0. If we drop the 
last and non-negative term in eq.(9), replace by in eq.(lO), and apply the 1-d WKB 
approximation, we obtain a lower bound of the classical action Sc- 

Slower = 2 r dzJ2m[VA{z) - Sq] , (11) 

with zi and Z2 the turning points which are the solutions of the equation Va{z) — So — 0- The 
tunneling rate calculated in this way gives the upper bound for the tunneling rate as pointed 
out in Refs. [10,11] This upper bound of tunneling rate corresponds to the physical situation 
in which the 2-d electrons follow the motion of the escaping electron instantaneously, the 
adiabatic limit, as a general theorem has shown.?? The classical trajectory corresponding 
the case in which r(T) = constant therefore 1 — cos[k • (r(T) — r(T')] = ( note B\\ex — 0). 
Setting exp{— A;2;(t')} = 1, integrating over r', and again using the 1-d WKB approximation, 
we obtain an upper bound for Sc from eq.(9): 

Supper = 2 dzJ2m[Vu{z) - So] , (12) 

with z[ and z'2 the turning points which are the solution of the equation Vu{z) — So = 0. 
Here the potential Vu{z) is 

/ X T. / N 1 /■^^ „ 27r/ce^no(l - A)^ , r , -> -.n9 
Vu(.)=V,U) + jI dk J^^^^ ' (exp{-fc} - 1)^ (13) 

This upper bound for Sc is smaller than that given by the frozen potential??, and is a better 
one. This is due to the partial inclusion of the dynamics of the plasma modes. 

The tunneling rate increases as Bj_^^ in the low field limit because of the suppression of 
the correlation potential discussed in the case of thermal activation rate. There is another 
reason for the increase of tunneling rate, the increase of the adiabaticity: If we fix the 
adiabatic potential Va{z), and allow the plasma frequencies to increase independently, then 
the e-e response will be able to follow the motion of escaping electron more closely, and 
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Sc decreases towards its lower bound, Sc — > Slower, the adiabatic limit. This result is the 
opposite of Ref.[10]. Because the damping is superohmic, the temperature dependence due 
to damping in the tunneling regime is found to behave as 0{T'^) and is small in the low 
temperature limit.?? 

Now we discuss the case that the parallel magnetic field is present, S||es 7^ and B^^^ — 0. 
By choosing the parallel magnetic field along the x-direction and = (0, — S||ej;2;, 0), 
according to the effective Hamiltonian eq.(5) we have the first two terms as 

(14) 



e 

P + -fill z 

c " . 



The x-direction motion is irrelevant. This Hamiltonian describes the 1-d motion of a particle 
coupling to a harmonic oscillator with a zero frequency. We integrate over the y coordinate 
to obtain the effective action for the motion in z direction as 

5e//[.(r)] =jdT ^-mz\T) + V^izir))] + ^ (^)' ^ J dr J dr' [.(r) - .(r')]^ . 

(15) 

This result suggests that the y direction motion acting as a dissipative environment to the 
z direction motion.?? The influence action from the plasma modes is ignored by assuming 
the adiabatic limit for simplicity, which will be discussed below. In the small field limit, 
the classical trajectory remains unchanged. The classical action can then be evaluated 
perturbatively: 
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(16) 



with semiclassical trajectory Zc{t) determined by the usual equation |mi;^ = Va{zc)—So. The 
tunneling rate decreases at zero temperature because of the bending of classical trajectory, 
and increases exponentially with temperature because of the excitation of the y-direction 
motion, with a fiyg^ dependence on the parallel magnetic field. In the calculation we have 
ignored a weak dependence of £0 on B\\ex- 
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We have discussed the crossing from thermal activation to quantum tunnehng, the effect 
of the correlation potential, the dynamical response of the 2-d electrons to the escaping 
electron, and the effect of weak parallel magnetic field. We now discuss the conditions 
to observe these effects. Let tol = \V^{zmax)\/fn be the small oscillation frequency in the 
inverse potential — Va{z). The crossing from thermal activation to quantum tunneling occurs 
at /cbTq = fi^b according to the standard theory.??. The frequency uji, also determines the 
dynamics of the escaping electron in the inverse potential. In order for the 2-d electrons to 
follow the motion of the escaping electron, the cutoff plasma frequency should satisfies the 
inequality ojc > ujb, and then the tunneling is essentially an adiabatic process??. Similarly, 
the condition for the weak parallel magnetic field B\\gx is eB\\ex/'>TT'C < To manifest the 
correlation effect in the escape process, we need ^yr^oZmax > 1 in the case of the thermal 
activation, and yJn^Z2 > 1 in the case of quantum tunneling, otherwise the correlation effect 
is small. 

Finally, let us make a brief comparison to experiments. The present theory is in agree- 
ment with experiments in the thermal activation regime, where consistent experimental data 
exist.?? The existing quantum tunneling data in the absence of magnetic field?? do not agree 
with each other and there is no satisfactory explanation. However, the magnetic field de- 
pendent tunneling results obtained in the present paper, the quadratic dependence on the 
parallel magnetic field and linear in temperature, agree with a recent experiment described 
in Ref.[5]. 
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